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o 

We consider a nonlocal version of the Nambu and Jona-Lasinio model. The non- 

(N ■ 

locality is contained in the quark-antiquark bilinears of the four-quark vertices as a 

a \ 

form factor of the Gaussian type. The model has three parameters which can be fixed 
in favour of the values of the pion mass, the pion decay constant /„-, and the current 
| | quark mass. The irir scattering amplitude is obtained by calculating the quark box 

and the cr-pole diagrams, where a is the scalar isoscalar meson. It is shown that this 

m ■ 

—i ■ amplitude satisfies the well-known Weinberg relation. We obtain the s,p,d wave 

m : 

scattering lengths in all isotopic channels and the s wave slope parameters. The 
results are in satisfactory agreement with both phenomenological data and the basic 
requirements of low-energy theorems, thus supporting to the form factor used. 

D 

43 ■ 



X 



1. INTRODUCTION 

The theoretical study of 7T7t scattering has a long-standing history. One of the pioneering 
works which has certainly boosted the development of chiral symmetric models was due to 
Weinberg [1]. In that paper the methods of current algebra [2, 3] were successfully used to 
compute the low-energy irn scattering lengths. After that many Lagrangian models of strong 
interactions of hadrons were intensively developed, using the linear [4] and nonlinear [5-7] 
realizations of chiral symmetry. Various low-energy properties of hadrons were described in 
the framework of these models, for instance tttc scattering (see e.g. [8-10]). 
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The chiral perturbation theory has been developed on the basis of nonlinear Lagrangians. 
This effective field theory describes the low-energy structure of different amplitudes in terms 
of expansion in powers of energies, momenta and current quark masses [11, 12]. Chiral 
symmetry determines the low energy behaviour of the 7T7T scattering amplitude to within very 
small uncertainties [13, 14]. For such a theory one takes for granted that chiral symmetry 
is spontaneously broken. There are theoretical arguments indicating that chromodynamics 
leads indeed to the formation of a quark condensate [15]. However, the chiral perturbation 
theory cannot explain this phenomenon. 

An analogy with superconductivity led Nambu and Jona-Lasinio (NJL) to the field the- 
oretical model of strong four-fermion interactions [16, 17] which provide for spontaneous 
chiral symmetry breaking and the formation of a quark condensate. The local version of 
this model was particularly well developed at the quark level [18-25]. The spectrum, low- 
energy dynamics, the main strong and electromagnetic decays, scattering and the internal 
characteristics of mesons have a reasonable explanation within this model. 

In this paper, we generalize the local approach to nn scattering by constructing the iyty 
scattering amplitude within the framework of a nonlocal NJL model with linear realization 
of chiral symmetry. The scalar and pseudoscalar degrees of freedom are considered. Let 
us note that one needs to know higher powers of the p 2 expansion to describe scattering 
lengths and slope parameters of ci-waves and higher. Such expansions were first a problem 
in the local approach: ad hoc form factors were inserted into the meson-quark vertices [26], 
their structure not being derived from the internal structure of the model, thus requiring 
justification [27]. 

A coherent method which takes into account the whole contribution of quark one-loop 
diagrams in the local NJL model has been formulated in [28, 29]. In particular, the method 
was successfully used to calculate nn scattering lengths. 

Recently, nonlocal versions of the NJL model have been actively developed [30-42]. These 
models are free from ultra-violet divergences, as well as there is hope to confine the quarks 
because the constituent quark mass depends on the momentum. The form factors of the 
meson-quark vertices are the direct consequence of this nonlocality. Thus, it seems inter- 
esting to calculate again nn scattering lengths and parameters of effective ranges in the 
framework of a nonlocal model. The relative simplicity of the process together with a large 
amount of experimental information about 7nr scattering serve us to gain a considerable 
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insight into the structure of nonlocal interactions from such calculations. 

The paper is organized as follows. In the next Sect, we consider the effective 577(2) x 
SU (2) symmetric nonlocal four-quark interaction. The dynamical mechanism of spontaneous 
chiral symmetry breaking is clarified after bosonization. We show that the model satisfies 
Goldberger - Treiman and Gell-Mann, Oaks, Renner relations. We fix parameters of the 
model in this section. In Sect. 3 we obtain the nn scattering amplitude, check Weinberg 
relation and calculate s, p, d scattering lengths in all isotopic channels and slope parameters 
of s-waves. The results are compared with empirical data. In the last Sect, we discuss them. 



2. NONLOCAL MODEL OF THE NJL TYPE 

The model is based on the effective SU(2) x SU(2) symmetric quark interaction in the 
scalar and pseudoscalar sectors 

C = q(x)(id x - m°)q(x) + | (J a (x)J a (x) + J^(x)J^(x)) , (1) 

where q, q are quark and antiquark fields, m are masses of current quarks, Tj are the Pauli 
matrices and G is the four-quark coupling constant. The nonlocal quark currents Ji(x) are 
given by 1 

Ji{x) = J d 4 x 1 d 4 x 2 5 (^x - Xl ^ X2 ^j f((xi - x 2 f) q(xi) T/ q(x 2 ), (2) 

where the function f(x) is normalized as /(0) = 1. The matrices Tj in eq.(l) are defined as 
follows r ff = 1, r£ = i^r a . 

One has for this Lagrangian after bosonization 

C = q(x)(id x - m c )q{x) ~ ^ {^ a {xf + a(x) 2 ) + J a (x)a(x) + 7r a (x)J^(x), (3) 

where a and n a are the scalar and pseudoscalar fields. The scalar field a has a nonzero 
vacuum expectation value (cr)o — °"o 7^ 0. To arrive at the physical scalar field with zero 
expectation value, it is necessary to shift the field: a = a + a . 
Varying the action with respect to a 

£>.-* 



1 The models with such nonlocality have been considered previously in [37-42] 
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one obtains the equation for the dynamical quark mass, the "gap" equation 

m(p 2 ) =m c + iG^f( P 2 ) J -0^J(k 2 )Tr [S(k)\ , (5) 

where Nf and N c are the numbers of quark flavours and colours, respectively. Note that 
due to the form factors the quark loop integral is finite. This equation can also be written 
as 

m(p 2 ) =m c - cr f(p 2 ) = m c + (m q - m c )f(p 2 ), (6) 

where m q is a dimensionfull parameter which plays the role of the constituent quark mass. 
As a result the quark mass becomes p 2 dependent. The quark propagator S(p) has the form 

p — m{p 2 ) 

The propagators of the meson fields are 

n ( *\ - 1 alAv 2 ) /ox 

aAP) -g-i + j„Ap 2 ) v 2 -mi; w 

where M CT)7r are meson masses, g^^ip 2 ) are the functions describing the renormalization of 
mesonic fields, and J CT7r (p 2 ) are the contributions of the quark one-loop diagrams 

J*AP 2 ) = J d 4 kf 2 (k 2 ) Tr [S(fc_)r CT>7r S(fc + )r CT)7r ] , (9) 

where k + = k+p/2, k^ = k — p/2. The meson masses M a>n are the poles of the propagators 

J.AK*) = ( 10 ) 

and the couplings g c<K {M 2 ^) are determined on the mass-shell of the mesonic states by the 
formulae 



dJaAP 2 ) 



dp 2 

The amplitude of the pion decay n — > \w is 



(11) 



A^(p)^ip"U, (12) 

where f w is the coupling of the weak pion decay 2 . The Goldberger - Treiman relation is held 
in the chiral limit f n = m q j 'g n (0) and /„- is given by the formula of Pagels and Stokar [45] 



oo 

N c f m(u) - ±um'(u) 



fn = T^2 / duum(u)—— — ■ (13) 



2 Detailed description of insertion of electroweak fields into nonlocal Lagrangian is given in appendix. 





Figure 1. Diagrams contributing to the amplitude of irn scattering. 

This expression is given in Euclidean domain u = —p 2 . 

The Gell-Mann - Oaks - Renner relation is fulfilled in the model 

M 2 f 2 = -2m c (qq)+0(m 2 c ), (14) 

where (qq) is the quark condensate. This relation can be obtained from the equation for 
the pion mass (10) with the help of the gap equation (5) and expanding in powers of M 2 
and m c . The form factor f(u) has been chosen in the form of a Gaussian exponent f(u) = 
exp(— u/A 2 ). 

The model has three parameters m c , A, G. After using the values of pion mass M n = 140 
MeV and the weak pion decay constant /„- = 92 MeV as input parameters we obtain instead 
of three arbitrary parameters only one m c . This parameter can be chosen in the interval 
4.05 — 7.425 MeV at the typical hadron scale 1 GeV [43]. Additional restriction can be 
obtained, e.g., from the ci-wave scattering length a 2 ,. As a result, the parameters of the 
model are m c = 4.8 MeV, A = 880 MeV and G = 22.87 GeV" 2 . Thus, one can find that 
m q = 345 MeV and {qq) = (— 255MeV) 3 . For the mass and the width of scalar meson we 
obtain M a = 570 MeV, IV = 220 MeV. 

3. 7T7T SCATTERING 



The amplitude of the process is described by two main diagrams plotted on Fig.l. 
Using the notation of Chew and Mandelstam one can write the matrix element of nn 
scattering 

— ^ii2^i3«4^(' s ? ti u ) + fiiwfiizuAit, s, u) + 5i x i A 5i 2 i 3 A(u, t, s) (15) 
in terms of three kinematical invariants 

s = (pi+p 2 ) 2 , t=(p 1 -p 3 ) 2 , u = (p 1 -p 4 ) 2 . (16) 
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Pi -/':.! ' 1 



For the amplitudes with the definite isospin values I = 0, 1, 2 one has 

T° = 3A(s, t, u) + A(t, s, «) + A(u, t, s) 

T 1 = A(t, s, u) — A(u, t, s) (17) 
T 2 = A(t,s,u) + A(u,t,s). 

The a-pole diagram contains the ann vertex which has the following structure: 

9^V\PIV\) = ^ / **/ ((* + |) 2 ) / ((* - f ) 2 ) / ((* 
xTr [S(A; + Pi)r CT S(A; - p 2 )T 7T S(k)T 7T ] , 

The box diagram contributes as 

A box (s, t, u) = 2^(M 2 ) (J 4 (s, t, u) + J 4 (s, u, t) - J 4 («, t, s)) 

j4(s , i ,„ ) . i »/^/((* + |) 2 )/((,-|) 2 )/((* + 

x / ^ + Pl ~ P 2 2 ~ P3 ^ j Tr [S(A;)r 7r S(A; - p 2 )r^S(A; + Pl - p 3 )r.S(A; + Pl )Y v \ , 

As a result, the whole 7T7T amplitude can be written as 

o 2 fs M 2 M 2 ) 

A(s,t,u) = A box (s, t, u) + Ag±{Ml) 9 °™\ Jj . (18) 

At the lowest order of the momentum and m c expansions one obtains the Weinberg result 

<? — M 2 

A(s,t,u) = —^ + 0(p i ) (19) 

J 7T 

Let us recall that in the local NJL model the Weinberg formula can be obtained on the 
basis of the box and the cr-pole diagrams, if one picks up only the divergent terms of the 
quark one-loop diagrams [26, 46] 



Pi-Ps' 2 



2 16m V s-M 2 
A(s, t, u) = - V + — -JL = _JL + ..., (20) 

1V1 IT b J IT 

where g = g w = g a , m is the constituent quark mass, and the sigma meson mass takes its 
model value M 2 = Ml + 4m 2 . 
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Scattering lengths af and slope parameters bf are the first terms in the expansion of 
partial-wave amplitudes for orbital angular momentum / and isospin I 

1 

in gatcm^ I dxP ^ x ) Tl ^ u ) = a i + k%I i + °^)- ( 21 ) 

-i 

Here s = 4(M^ + k 2 ), t = — 2/c 2 (l — x), u = —2k 2 (l + x); k is the cm. three-momentum, 
and x = cos 9, where 9 is a scattering angle between pions in the cm. frame and Pi(x) are 
the Legendre polynomials. 

To obtain scattering lengths we expand the integrand of eq.(21) in powers of k 2 . The 
subsequent integration over x lead to the expressions for the scattering lengths and slope 
parameters that are given in appendix. 

The obtained values for the scattering lengths of the s,p and d waves in all isotopic 
channels as well as the parameters of the effective ranges for the s waves are given in Table. To 
compare our results, we also present there characteristics obtained in the framework of the 
local NJL model [28]. The empirical data are taken from the papers [49] (the life-time of the 
7r + 7r~ mesonic atom), [50] (the K — > 3n decay), [48] (the K ± — > -K + it~e ± v e (u e ) decay), and 
[47] (the phenomenological analysis of data). 

To check the result, we passed to the local limit in our formulae, by regularizing the 
quark one-loop diagrams in the standard way and taking A — > oo in the nonlocal form 
factors. After this the results of known calculations in the local model are reproduced [28]. 



4. CONCLUSIONS 

The present calculations show that the main low energy theorems of pion physics are ful- 
filled within the nonlocal NJL model extension considered. We have checked the Goldberger 
- Treiman relation, the Gell-Mann, Oaks, Renner formula and the famous Weinberg result 
for the low energy tttt scattering amplitude. 

The Goldberger - Treiman relation is an exact result of the theory in the Goldstone limit 
and reflects the chiral asymmetry of the ground state. In the NJL model as well as in the 
linear sigma model, the Goldberger - Treiman relation can be understood as (<Jo) = f n , 
where (a ) represents the point at which the effective potential has a minimum. 

The Gell-Mann, Oaks, Renner formula in the NJL model relates the quark condensate 
with the gap equation, thus reflecting the dynamical mechanism for the spontaneous chiral 
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Table 1. Scattering lengths and slope parameters for the irir system. Colums LNJL, PY and EXP 
correspond to the results obtained in the local NJL model [28], results of the phenomenological analysis of 
experimental data [47], and known data extracted from the life-time of tt + it~ mesonic atom [49], from the 
K -> 3tt decay [50] and if ± -» 7r+7r"e ± ^ e (^ e ) data [48]. 





This work 


LNJL 


PY 


EXP 




0.20 


0.190 


0.230 ±0.015 


0.216 ±0.017[48] 


alM n 


-0.049 


-0.044 


-0.048 ± 0.0046 


-0.0454 ± 0.0049 [48] 
-0.041 ± 0.036[50] 


(ag - a 2 )M n 


0.25 


0.234 


0.277 ±0.014 


0.264™ [49] 
0.268 ± 0.027[50] 


a\ x 10~ 3 M3 


37 


34 


38.4 ±0.8 




ajj x 10~ 4 M5 


16.3 


16.7 


18.7 ±0.41 




a\ x 10~ 4 M5 


2.8 


3.2 


2.78 ±0.37 




blMl 


0.26 


0.27 


0.312 ±0.014 




blMl 


-0.085 


-0.079 


-0.090 ± 0.006 





symmetry breaking. 

The nn scattering amplitude has been obtained in the Hartree-Fock approximation of 
the nonlocal NJL model and was used to calculate the threshold characteristics. The results 
for the scattering lengths in s,p and d waves and for the s-wave slope parameters are in 
satisfactory agreement with the empirical data. 

The main contribution to the values of the s,p scattering lengths is given by the p 2 
terms of the amplitudes. By virtue of the Weinberg relation, they are less sensitive to the 
structure of the form factor. The d- waves are mainly determined by the p A order terms of 
the amplitude and the sufficient predictions for the ci-waves are a non-trivial result of our 
calculations. In terms of box and cr-pole diagrams, they are additive for the a® scattering 
length and have opposite signs in the a\ case. The value of the a\ scattering lengths is 
strongly affected by the magnitude of the a meson mass predicted by the model. 

Returning to the problem of fixing model parameters, note that the best agreement with 
experimental a\ is achieved for m c = 4.8 MeV(a2 = 2.8 x 10 _4 M~ 5 ). For the current quark 
mass equal to 4.4 and 5.2 MeV this scattering length equals 3.5 x 10~ 4 M" 5 and 2 x 10 _4 M- 5 , 
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respectively. 

In the future we plan to extend the nonlocal model considered by including the vector and 
axial- vector mesons and taking into account the 7r — cli transitions. It also will be interesting 
to consider the U(3) x U(3) chiral symmetry group and study the internal structure and 
interactions of the pseudoscalar, scalar, vector, and axial-vector nonets. Of special interest 
will be processes in which the p 2 dependence of the amplitudes will be relevant. 

Moreover, it is interesting to investigate the nn scattering process in the hot and dense 
medium in the framework of a nonlocal model. When chiral symmetry is restored, the a 
meson mass tends to the mass of pion [42, 46, 51]. This can lead to a noticeable increase 
in the irn scattering amplitude near the phase transition of hadron matter into quark-gluon 
plasma. Such investigations in the framework of local models are performed in [52]. 
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APPENDIX 



Weak pion decay constant 

To determine the weak pion decay constant one inserts external sources associated with 
the electroweak fields. For this purpose it is necessary to take into account the contributions 
induced by the nonlocal interaction. We use the standard procedure: the external fields are 
introduced via the Schwinger phase factor which delocalizes quark fields 

q(y) -> Q(x,y) = Pexp ji j dz»{V«(z) + ^) 75 ]T a j q(y). (22) 

Here V is the Dyson ordering along the path, V£(z) and A^(z) are the external vector and 
axial-vector gauge fields which can be identified with electroweak fields, T a stands for the 
generators of the flavour group. Since the derivative of the path integral of any function 
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k 




k- 



Figure 2. Diagrams corresponding to the weak decay of the pion 7r — ► fiv. 

Fp(z) over the path does not depend on the trajectory along which one calculates the integral 
[44] 




As a result, non-minimal interactions, which could be induced by the kinetic term of quark 
fields, are absent. The gauge invariant procedure used here leads to additional vertices 
with any number of external fields which have their origin in the nonlocal nature of the 
interaction. 

Weak pion decay is described by three diagrams (see fig. 2). Note that in the chiral limit 
the contributions of nonlocal diagrams are exactly cancelled. 



Scattering lengths and slope parameters 



a2 ° = 3^(^(MM^0) + A(0,0,4M2; 
(«g-«S) = ^1(4^,0,0) 

a{ = ^^(V)(0,4M^0) + A( S ^ 

-2A (a,t) (0, 0, 4M*) + A (s ' s) (0, 0, 4M^) + A (s ' s) (0, 4M^, 0) 
(a° 2 -a 2 2 ) = ^^(^»>(4iU^0)-2^ 
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bl = Y^r( 2A(U) (°' ' 4A^) + 2A(*)(0,4M2,0)-AW(0,4A^,0) 



(0, 0, 4M^) - A w (0, 4M^, 0) - A {s) (0, 0, 4M^) 
(6g — 6g) = Y ^(2A( s \AM^,0,0)-A( t \AMl0,0)-A^\AM 2 n A0)y 



where 



O^ 1 
A {x \a,b,c) = —A{s,t,u)\ s=a ,t= b , u =c, A (x ' y) {a,b,c) = g^ A ( s , *, u)\*=a,t=b,u= 
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